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1.0 INTRODUCTION 

The data from the Sky lab Earth Resources Experimental Package 
(EREP) is helping to answer questions now posed by many remote sensing 
specialists and engineers regarding the sensor specifications of 
future earth resource satellites, for example, using mul tispectral 
scanner data from the Earth Resources Technology Satellite (ERTS-l), 
we can now identify many important agricultural crops with a high 
degree of accuracy. Therefore, the two obvious questions that arise are 
(1) Will we be able to identify agricultural crops any more accurately 
with the EREP S 1 92 mul t i spect ral scanner, which nominally has three 
times the resolving capabilities of the scanner on ERTS-l? and, (2) 

Will this higher resolution enable us to reduce the number of sub- 
samples required to achieve a particular degree of crop identification 
accuracy? 

Another aspect of the research to be done with Sky lab EREP data is 
that of technique development. Due to the complex nature of the sensor 
(i.e., the 13 bands from the S 1 92 mul t i spectral scanner and the seven 
photographic products from the S190A and S1908 camera systems), much 
research must be done to determine what data are the most relevant to 
monitoring agricultural resources and how that data can be handled 
most efficiently. 

The emphasis of the Skylab agricultural investigations being car- 
ried out by the Center for Remote Sensing Research (CRSR) is placed on 
evaluating quantitatively the Skylab data with respect to its useful- 
ness in agricultural inventories. In addition, we will attempt to 
determine how well these results compare with those that can be 
obtained from ERTS-l data. The following sections of this quarterly 
progress report explain in detail (1) what has-been accomplished to 
date and (2) those tasks which are projected for the remainder of the 
contract period. 
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2.0 WORK PERFORMED DURING THE PERIOD COVERED BY THIS REPORT 

2 . 1 Status of Skylab 2 and Supporting High Altitude Aircraft Imagery 

The only Skylab 2 overpass covering our agricultural test sites 
was Track 63 made on 'June 3, 1973* At this time, we have received 
all photographic products from the S190A mul t ispectral camera system 
and the screening images for the S 1 92 mul t i spectra 1 scanner. The 
corrected S 1 92 tapes have been ordered, but as of yet have not been 
received. No S 1 9 0 B imagery was obtained of the test sites. 

Aircraft photo missions of the agricultural test sites have been 
flown to coincide as closely as possible to the overpass of Skylab 2. 

To date, high quality RC-8 color photographs and 70 mm multiband 
imagery have been recei ved from the Santa C 1 ara Val 1 ey , the Salinas 
Valley, and the west side of the southern San Joaquin Valley test 
sites. 

2 . 2 Agricultural Test Sites 

Due to the westerly shift of the nominal orbit of Skylab 2, only 
partial coverage of the two primary test sites was obtained. For the 
San Joaquin County test site (#101-106), the southwest corner was 
imaged by the S190A mul t ispectral camera only. For the southern San 
Joaquin County test site (#101-104), the western half was imaged by 
the S190A and a narrow strip along the western border was imaged by the 

si 92. 


In view of this limited coverage of the primary test sites, the 
existing sites were expanded, with verbal agreement by personnel at 
Johnson Spacecraft Center (JSC) , to include the agricultural areas in 
the Santa Clara Valley and the Salinas Valley (see Figure 1). These 
areas will also be imaged by the Skylab 3 mission thereby providing 
the sequential coverage that is needed for many of the proposed inter- 
pretation tasks. 

2.3 Ground Data Collection 

Thus far,. most of the efforts of the CRSR personnel have been 
concentrated on collecting ground data in support of the Skylab agri- 
cultural experiment. These data are required for the training and 
evaluation of both human and automatic image analysis. Unfortunately, 
the data that were collected in the two primary test sites -- San 
Joaquin County and southern San Joaquin Valley -- prior to the launch 
of Skylab 1 will be of little use in this experiment due to the 
unplanned shift in the Skylab 2 ground tracks. However, we were given 
sufficient notice of this shift to obtain low level aerial photos and 
to collect ground data in the Santa Clara Valley and the Salinas Valley 
nearly coincident with the June 2, 1973 Skylab overpass on Track 63. 
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Figure 1. The iocatjons of the new agricultural test sites — Sa 
Clara Valley and Salinds Valley — with respect to Skylab ground 
tracks. (Solid lines delimit test s.Jtes and dashed lines delimit 
S190A photographic coverage.) 



2.4 Personnel Training 


An internal training session is being given at the CRSR to image 
analysts to familiarize them with the CALSCAN* automatic classifier. 
Prior to this time and when working with ERTS-1 digital tapes, the 
task of image stratification, which is a prerequisite for efficient 
computer analysis, was done by human image analysts. The stratified . ^ 

imagery was then given to a computer operator who trained and classified 
the individual strata' with data from the digital tapes. It was felt, 
however, that these two tasks could be done most efficiently by the 
image analyst alone because, with his greater knowledge of agricultural 
cropping practices, he can more effectively train the CALSCAN classi 
fier and make a better evaluation of the final results. 

2 . 5 Automatic Interpretation 

The planned automatic data processing studies call for extracting 
agricultural information from EREP computer-compatible tapes and pho- 
tographic transparencies. The processing of both types of data will 
give a comparison of classification results using an analog storage 
medium (film) versus digital storage (tape). 

To date we have not received the digital tapes, but processing of 
the photographic transparencies has begun. The northern half of the 
Salinas Valley test site Is being used for preliminary evaluation of 
the technique. This area has been scanned with a microdensitometer on 
all four bands of S190A photography. The scan interval along the X and 
Y axes was .01 inches, and the aperture of the scanner was set at .001 
inch diameter so that each data point represented a 1.12 acre spot on 
the ground. The resulting density measurements were recorded on mag- 
netic tape. This tape is now being reformated on the CDC 6600/7600 
system so that it will be compatible with the CALSCAN program. 

3.0 WORK PLANNED DURING THE NEXT REPORTING PERIOD 

3 . 1 introduction 

Due to the uncertainties that have been associated with the EREP 
missions, it has been difficult to specify those tasks which will be 
accomplished throughout the remainder of the contract period. However, 
now that the areas and time of Skylab coverage are known, we can pro- 
ject fairly accurately those tasks which will be done and their appro- 
priate time frames. A summary of these tasks v is shown in Figure 2. 


* CAL SCAN is the CRSR version of the LARS-Purdue pattern recognition 
program adapted to the CDC 6600/7600 system at the University of 
California, Berkeley. 
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2 . Training 

3. Strat i f i cat ion 
k. Irrigated Study 

5 . Automat i c Interpretation 

a. photos 

b. tapes 

6. Manual Crop Inventory 

7. Multistage Sample 


















3.2 Agricultural Land Use Classification 

The stratification of agricultural land from other land uses is 
the first step for all analyses of Skylab data with respect to agri- 
cultural resources. Such information is not only necessary for mon- 
itoring land use change and the efficient allocation of ground samples, 
but also is a vital input to the CALSCAN classifier so that resulting 
statistics can be obtained in the most cost effective manner, without 
affecting classification accuracy. 

i 

The stratification of agricultural lands will be done in all test 
sites. The bases of these stratifications will be on (1) visual appear- 
ance, and (2) the classification system that is now used by the Calif- 
ornia Crop and Livestock Reporting Service (CLRS). The first strati- 
fication system which will be input into the data bank for future use 
by CALSCAN, will be more detailed than the second. The second strati- 
fication is the type that would be used to allocate ground samples for 
a conventional enumerative survey. 

The accuracy of the visual stratification system will be deter 
mined indirectly using the CALSCAN classifier. Several different 
versions of the strat i f icat ion boundaries, e.g., boundaries drawn by 
different image analysts, on different dates, on different film types, 
and on ERTS-1 composites, will be input into CALSCAN to determine how 
they affect the accuracy of the final output. 

The second stratification system will be evaluated by comparing 
it with the existing strat i f i cat ion boundaries as delineated by the 
CLRS using HAP IT. 

3 . 3 I rrigated Land Study 

There presently exists an important need on the part of the 
Department of Water Resources of the State of California for a periodic 
tabulation in any given year of the statewide acreage of agricultural 
land receiving irrigation. A study will be performed to investigate 
the extent to which this tabulation can be accomplished using two dates 
of EREP data and the appropriate sampling designs. This study will not 
be able to give a valid yearly statistic of the irrigated lands within 
the test area due to the fact that only two dates of Imagery area 
available and imagery taken at other dates would be required. How- 
ever, it will help determine (l) the feasibility of making such 
inventories, (2) the expected accuracies of such inventories, (3) if 
similar techniques could be used on ERTS-1 imagery, and (4) a projected 
cost for making statewide inventories. 
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3,4 Automatic Interpretat ion 


3.4.1 Photographic Data 

After the reformatted tapes are available (see section 2.5), the 
data obtained from- the scanned S190A photos will be analyzed with the 
CALSCAN classifier for purposes of crop identification. If the resul s 
of this analysis indicate that the results are comparable to those 
expected from the SI 92 mul ti spectra 1 scanner, the remaining agricultur 
areas not imaged by the scanner but by the photos will be scanne wi 
the microdensitometer. 


3 . 4.2 hul t i spect ra 1 Scanner Data 


The initial experimentation with data from the S 1 92 multispectra 
scanner will be to determine which combinations of the 13 channels are 
optimum for identifying crop types. This analysis will be made wit 
the CALSCAN classifier. 


3.5 Manual Crop Inventory 

To date most of the identification aiid inventory of agricultural 
crops using satellite data (i.e., ERTS-l) have been done semi- . 
automatically by computer analysis of digital tapes with human inputs 
of training materials. Although such systems have achieved very accu- 
rate results in certain test areas, there is still a need to develop 
techniques for manual interpretation of satellite data for agricultural 
resource inventories. While computer-based systems may ultimately 
prove to provide the most efficient method for gathering agricultural 
statistics of extensive areas, at the present time the human inter- 
preter represents the most expedient way to perform an operationa 
inventory in the United States. For many of the emerging nations ot 
the world where both national agricultural statistics and computer 
systems may be non-existent, the data gathered by human interpreters 
can provide a valuable input to the management decisions for agricultural 

resources. 

In any attempt to develop efficient techniques for the human 
interpretation of Skylab data, several factors must be considered: (1) 

because of the large areal coverage of the Imagery, 100 percent image 
interpretation of the entire frame for detailed Information is not 
practical, (2) a simple method is needed to evaluate the accuracy of 
the interpreter ' s estimates and, if necessary, to adjust these esti- 
mates, and (3) the low resolution of the imagery makes accurate acre- 
age estimates by human interpreters impossible. 

An inventory technique employing a double sampling design util- 
izing point data is expected to deal effectively with these constraints. 
For the first stage (large sample) the interpreter will determine the 
presence or absence of the crop of interest at each of a large number 
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of points throughout the survey area. These data will then be used to 
estimate the proportion of the area that is planted to that particular 
crop. This proportion, when multiplied by the total area being 
inventoried will give an estimate of total crop acreage. The second 
stage in the inventory will consist of a subsample of the large samp e. 
These subsample points will be field checked, and the correlation 
between ground conditions and image interpretation estimates w» 1 1 be 
used to evaluate interpreter accuracy and to calculate a ratio esti 
mator to adjust the interpreter's estimated proportion. . It must be 
noted that, due to the long delay between data acquisition and the 
receipt of EREP imagery, the second stage data will probably have to 
be collected from high altitude photography that was obtained coinci- 
dent with the Skylab overflight. 

3.6 Multistage Sampling of Agricultural Resources 

In general, there are two types of basic information required by 
agricultural planners: (l) an estimated the quantity of specific 

types of agricultural resources in each administrative unit, and (2) 
an inplace map of those resources. When a high correlation exists 
between estimates of the resources as made from spacecraft data and 
those based on ground observations, unique and va 1 uabl e i nforma t i on 
for meeting both of these requirements can be provided in a cost- 
effective manner. 

The culmination of our Skylab agricultural experiments will be 
the demonstration of agricultural survey techniques using a multistage 
sampling model. Through use of the d i scr i mi nant ana 1 ys I s techniques 
described below, Skylab photo and scanner data will be combined with 
aerial photos and ground data to give an estimate of crop acreage 
within the areas common to Track 63 of Skylab 2 and Track 59 of Skylab 
3. We have not yet chosen which crop or crops will be inventoried; 
this will be determined after the data from the S 1 92 mul t i spectral 
scanner have been analyzed. Under normal operational conditions, the 
optimum time for obtaining remote sensing data for a successful agri- 
cultural inventory is known and is carefully adhered to. However,, 
since this is not possible under the restrictions of the EREP missions, 
the crops to be surveyed wil 1 be those which can best be inventoried 
with the available Skylab data. 

The multistage model relies heavily upon the first stage in which 
the information extracted from the Skylab data. by human interpreters 
and automatic classifiers provides the ini tial estimates of the . resource . 
The first step of the data extraction process will be one in which 
human interpreters stratify all fields within the area of interest 
into broad land use categories and crop classes based on the! r . appear- 
ance on the Skylab SigOA ektachrome imagery. At- this time political 
and geographic boundaries will be superimposed on the imagery to fur 
ther define the areas of interest. Next, a number of fields which 
represent the various agricultural resources of interest will be 
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selected from each stratum to train the discriminant analysis pro- 
gram. The identities of these fields will be determined from ground 
data and/or the interpretation of aerial photos. The number of train- 
ing fields required for each crop class will be dependent upon the 
variability of the spectral signatures of the crops present. This 
variability is caused by such factors as different cropping practices, 
local soil differences, and genetic variations within each particular 
crop type. For example, a crop such as alfalfa which may be in severa 
stages of maturity throughout the survey area at the time of image 
acquisition may require five or more fields per stratum for adequate 
training, whereas only one training field per stratum may be needed for 
a less complex crop type such as corn. After the fields have been 
chosen, they will be located on, and extracted from, the spacecraft 
imagery. The mu 1 1 i spect ra 1 data from the training fields will be run 
through the discriminant analysis program to obtain a po i nt -by-po i nt 
classification of the entire area by strata. This procedure will pro- 
vide an initial estimate of the acreage of the vegetation classes. 

In second stage of the model, the results of the discriminant 
analysis will be sampled to determine their relationship to ground 
estimates of the resource. Sampling units (SU) will be defined by 
dividing the entire area into rectangular blocks. The size and shape 
of these blocks will be determined by (1) the information requirements, 
(2) the variability that occurs in the SU estimates as the block size 
changes, ( 3 ) the costs of further subsampling, and ( 4 ) the resolution 
of the scanner imagery. 


Probability sampling is expected to provide the most efficient 
sampl i ng design that can be appl ied in the second stage of the model . 
Probability sampling is a special case of the mean of the ratios esti- 
mation where samples are allocated proportional to the expected va r th 
iance of the X. estimate. For this model, the total value of the t 
SU, denoted by 1 X . , is evaluated by 


X. 

i 


E £ 

m=l j=l 


I V. 
m j 


where 



M 

V . 
J 

J 


lifC = j 
m 

0 otherwise 

Am L 

crop class for the n "pixel*' (picture element) 
of the SU, as determined by the discriminant 
anal ys i s , 

the number of "pixels" per SU, 
the crop class being evaluated 

the number of crop classes. 
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The value or weight (v.) is assigned to rank the various crops or veg- 
etation types based on J their relative importance to the survey. In an 
agricultural inventory where total dollar value is the objective, the 
v.'s are the average dollar values per “pixel 11 of the various crops 

(j). tf the making of an inventory of a single crop is the objective, 
the value of the crop of interest is 1 and all other v.'s are set 
zero. In many cases this weighting factor is primarily affected by the 
marketing conditions of each crop, and is highest for those crops for 
which acreage estimate errors are most important. 

The variance of the population is estimated by 

s2 = ^X>i -^ 2 


The number (n) of SU ' s to be selected for photo and ground measure- 
ment when no remote sensing information is available is determined by 


n 


Nt 2 S 2 

N (AE) 2 + t 2 $ 2 


where AE - the allowable error, in units of value 

t is a value obtained from "students t" tables and 
2 

S is as defined previously. 

The n points are then selected from the list of SU 1 s proportional 
to their estimaged value. 

The selected SU's are then carefully transferred to the correspond- 
ing high flight photography where precise field size measurements are 
made for use later in adjusting the acreage estimates obtained from the 
class? f ier. 

From high flight images, low altitudes images, ground identifica- 
tion and historical data, the "correct" classification for each field 
in the SU is determined, down to crop type and maturity. 

The total value for the area (T) is estimated using the probability 
of selection (P.) and the photo/ground estimate of SU value (Y.) by 
means of the relation, 
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, n 

T.iS 

n i = l 


Y. 

i 


P. 

i 


x. 


where P. - 


N 

ST' X . 

i=i 1 


The variance of the estimate for T is 


s „ 2 = Var(T) = “ ^ p ; ^ 
T n i=l i 


T) 


2 


If the photo/ground estimate (Y. ) were to be perfectly proportional 
to the remote sensing estimate (X.), only one ground sample would be 

needed to determine the proportionality constant. More realistically 
however, the number of ground samples (n) for future surveys is esti 
mated by: 


2 2 2 
N (AE) + t Sy 


This probability sampling model is appropriate when a single 
parameter such as “acreage of a single crop", "value of all the crops 
present", or "demand for irrigation water" is desired. Such a model 
can be replaced by a regression sampling model if estimates on a cropr 
by-crop basis are required; however, the regression sampling mode will 
only meet the allowable error criterion for the total value of all 
crops by strata. 



